Radiative and leptonic decays of the pseudoscalar charmonium state rjc 
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The radiative and leptonic decays of rjc 77 and rjc — » l^l^ are studied. For 77c — > 77 decay, 
tiie second-order electromagnetic tree-level diagram gives the leading contribution. The decay rate 
of rjc 77 is calculated, the prediction is in good agreement with the experimental data. For 77c — > 
I'^l^ , both the tree and loop diagrams are calculated. The analysis shows that the loop contribution 
dominates, the contribution of tree diagram with intermediate state can only modifies the decay 
rate by less than 1%. The prediction of the branching ratios of rjc — > e+e" and i-i^f^" are very tiny 
within the standard model. The smallness of these predictions within the standard model makes the 
leptonic decays of r/c sensitive to physics beyond the standard model. Measurement of the leptonic 
decay may give information of new physics. 

PACS numbers: 13.20.Jf 



The radiative and leptonic decays of rjc "f"f and 
rjc — > involve electromagnetic and weak interactions. 
The decay rates are determined not only by electroweak 
interaction, but also by strong interaction, which binds 
the quark-antiquark pair cc in rjc together. The decay 
amplitudes of these decay processes are generally convo- 
lutions of the wave function of rjc and the electroweak 
transition amplitude of cc 77 and/or cc . The 

decay process of rjc 77 can be used to test the decay 
constant and the wave function of cc in this meson state. 
The decay rate of 77c —>■ 77 has been measured in experi- 
ment Within the framework of the standard model, 
the leptonic decay rjc — > is contributed by fourth- 
order electromagnetic transition and tree-level weak tran- 
sition induced by Z^. The low probability of a fourth- 
order electromagnetic and weak transition makes the lep- 
tonic decays sensitive to hypothetical interactions arising 
from physics beyond the standard model, such as the ex- 
istence of a light pseudoscalar Higgs boson in the next- 
to- minimal supersymmetric standard model [1, 0, 0, S 
or leptoquark bosons that carry both quark and lepton 
flavors [3, H, S J both of which can enhance the lep- 
tonic decays of rye with appropriate values of new physics 
parameters in these models. 

In the literature there have been many theoretical 
works on the calculations of two-photon decays of pseu- 
doscalar heavy quarkonia rjc, rjh etc., based on relativis- 
tic quark model or potential model [ill, [13, [I3l , Bcthe- 
Salp eter equation [ijj ITal. heavy-quark spin symmetry 
[la |. and Lattice QCD [l3|. However the decay rates of 
leptonic decays of pseudoscalar heavy quarkonia have not 
been known yet. 

In this work the radiative and leptonic decays of 77c — > 
77 and r7c are studied consistently within the 

framework of the standard model, by using a method 
different from those used in the literature. The effec- 



tive Hamiltonians in quark level for these decays are cal- 
culated at first. Then the factorization formula is de- 
rived, the decay amplitudes are expressed as convolutions 
of meson wave function and the hard transition ampli- 
tudes, where the wave function is controlled by nonpcr- 
turbative QCD, for which I use the result calculated in 
QCD sum rule [l8| . while the hard transition amplitudes 
are calculated with perturbation theory. For the process 
?7c l^l^ , the loop diagrams are calculated analytically. 
The infrared divergence is analyzed, the possible infor- 
mation of new physics is also briefly discussed. 
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FIG. 1: Diagrams for 77c 77. 



The diagram for the decay of 77c 77 is depicted in 
FigU The transition matrix element relevant to Figl^a) 
can be written in coordinate space as 



(2^ 



(1) 
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where e is the absolute value of the charge of electron, Qc 
the charge of the c quark in unit of e, A''(a;) and A'^{y) 
are the electromagnetic fields for the two photons. Con- 
tracting the creation operators of the two photons with 
that in the electromagnetic fields, then the amplitude Ti 
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becomes 



(2^ 



-ik-{x-y) 



*tJ- ikix *i/ ik2y 



(2) 



/c — m, 



-l^ceiv\al3^i e " £2 e 



•(0|c(a;)ac(y)/3|77c), 

where a and /3 are the Dirac spinor indices. It is un- 
derstood that the repeated indices are summed. In both 
eqs. ([ij and ([2]), summations over color degrees of free- 
dom and a factor Pexp{i / dz'^T°'A'^{z)} between c{x) 
and c(j/) that makes the amphtude gauge-invariant are 
indicated. 

The leading-twist wave function for rjc meson can be 
defined through the matrix element {0\c{x)ac{y)p\r]c) [H, 
[20l |. which will be further discussed in Appendix A, 



{0\c{x)ac{y)f3\r]c) 



(3) 



where u = 1 — u, /i is an energy scale, and frj^ is the decay 
constant of 77^, which is defined as 



(4) 



here is the four-momentum of the meson. 

With the matrix element given in eg.®, the transition 
matrix element Ti becomes 

Ti = -\uQy I' du J d^xdS J ^ (5) 



Then the total contribution of Fig[l] (a) and (b) is 

A^Ai + A2. (9) 

Using fc^ = fc| = and 2ki ■ k2 = , one can simplify 
the amplitude A to be 



A = 2if„^Qy f du 
Jo 



uumi + Toi 

Ijc ^ 



(10) 



For simplicity in the following calculation, one can define 
a new quantity M as 



M = 2ifr,^Qy I du 



(t){u,ii) 



then we can obtain the square of the total amplitude 



(11) 
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\Af = ^ml\M\'' 



(12) 



The decay width can be calculated by the following for- 
mula 



1 1 I J lO 1^1 

2!8^' ■ mt 



(13) 



where ^ is the statistic factor for two identical particles, 
k is the three-momentum of one of the photons in the 
rest frame of rjc- 

Next let us discuss the leptonic decay j]c ^ l^l - In 
the standard model, the Feynman diagrams for this pro- 
cess include photon-Fermion loop diagrams and tree 
diagram, which are depicted in Fig[21 



It is not difficult to perform the integration over the co- 
ordinates X, y and the momentum k. After these manip- 
ulations, one can obtain 



Jo 



du- 



cj){u,n) 



(uki — uk2)^ ~ 
■ep,^,kPk^4''en2nfS\k, + k^-p). 



(6) 



Remove the overall four-momentum conservation factor 
(27r)"'(5^(fci + k2—p), one can obtain the decay amplitude 
contributed by Fig[T]^a) 



Jo 



du 



(t){u,fi) 



{uki — uk2Y — ml 



(7) 



The contribution of FiglTJb) can be obtained by making 
the exchange of fci ^2 and ei -f-> £2 in Ai 



A2 ^ Ai 



ei ^ €2 
ki ^ k2 



(8) 




FIG. 2; Diagrams for rjc — > (a) and (b) QED contri- 

bution where two virtual photons as intermediate state, (c) 
weak Z° contribution. 



The effective Hamiltonian describing the transition of 
cc l^l^ can be calculated with the help of the Feyn- 
man diagram in Fig[21 The contribution of Fig[2] (a) is 
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calculated to be 



2„4 



Hi = Qie 



(2^ 



(14) 



{(A;2 - mf + ie)[{k ~ k^f + ie][(fc + ^2)^ + ie] 



x[(fc-/ci+pc)^-TOc+*e]}, 

where mi is the lepton mass, pc the four-momentum of 
c quark. The quark fields c and c are in the momen- 
tum space. They are related to the field operators in the 
coordinate space by 



(15) 
(16) 



Then the effective Hamiltonian }i\ can be reexpressed in 
terms of quark fields in the coordinate space 



c= / d4a:c(a;)e*P^ % 



Hi=J d'xe^P-^J d'ye^P-yQyJ ^ 



(17) 



where au is defined as 



ak = (k^ -mf + ie)[{k - kiY -f ie] 

X [{k + k2f + ie] [{k - ki + pcf - ml + ie 



(18) 



With the effective Hamiltonian Hi, the amplitude con- 
tributed by Figl2] (a) is 



(27r)4 (27r)4 



(19) 



(2^(2^ 



d'^k 



uiki)-f^i/: + mi)Yvik2) 



The matrix element in the above equation can be treated 
as 

{0\c{x)j^{/:- i>c + mc)7,,c(y)|?7,) (20) 

= {Q\c{x)ac{y)p\iic)Yfy{Ji- Jki+ i>c + mc)^^]al3. 

With the help of eq. and the wave function of rjc de- 
fined in eq.Q, the decay amplitude Ai can be expressed 
as 



2„4 



d'^k 
(2^ 



u{ki)Y{^ + mi)Yv{k2) (21) 
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Perform the trace operation and using the identity of the 
gamma matrices 

7q7/37A = OaPlX + 9l3\la " 9aXlp + ie^a^x^^^^, (22) 
the amplitude Ai can be further calculated to be 

Ai^tQ^f^f du(t>{u,fi)[u{ki)rj5 (23) 
Jq 

■ v{k2)aia - miu{ki)a'^'-''v{k2)e^iypap''a2], 



where 



a-ic 



d^k 



(24) 



Pak • (fc — fci + up) — (fc — fci + Up)aP ■ k 



a-k 



d^k {k-ki+upY 
(2^ '^k 



(25) 



ak 



Next one needs to perform the loop integrations in the 
coefficients aia and 02- 

Using p ~ ki + k2 and a few steps of algebra manip- 
ulation, one can cancel one propagator in the coefficient 
aia, then the four-point loop integration of aia can be 
reduced to a sum of a few three-point loop integrals 



d^k 
(2^ 



{ Pa 



1 

[(fc - fci)2 + ie] [(fc + fc2)^ + ie] [(fc - fci + upY - ml + ie] 
u 



(26) 



(fc2 — m^ + i^)[{k — kiY + iA[{k — fci + upY — ml+ ie] 
(fc — uki + uk2)a 

(fc^ — mf -f ie)[{k — kiY + ie][{k — fci + upY — rnf. + ie] (fc^ — mf -f ie)[{k + fc2)^ + ie][(^ ~ fci + upY — in-f. + ie] 

I 



(fc2 - mf + ie)[{k + k2Y + ie][{k - ki + upf - ml + ; 

(fc — uki + Uk2)a 



With the Feynman parameterization, the integration over the loop momentum in the above equation can be 
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performed, then the coefficient aia can be expressed as 
integrals over Feynman parameters 



aio 







1 pl-x 

dx I dy{pa 



1 u 



(27) 



u 
as 



[y + u{x — l)]fciQ + uxh 



2a 



as 



-uxkia + (u — ux — y)k 



la 



02 



with 



fli = rn^^lux + y){ux + y — 1) + xuu] + xm^ (28) 

- 

02 = [\ — 2x — 2y + [x + ?/)^]to^ + xu[xu — y)Tn^^ 

- 

03 = [1 — 2x — 2y + {x + y)'^]mf + xu[xu — y)rn^^ 



While the four-point loop integration of the coefficient 
in eq. (|25p can not be reduced, it can be expressed as 
a 3-fold integral of Feynman parameters 



167r2 



1 pl—x pl—x—y 

dx I dy I dz 











(29) 







where 



{ux + z — u)k1 — (ux + y ~ u)k2 



(1 — X — y — z)'^'m^ + [xuu 

— lux + z){ux + y)]m^^ + xm^ — it. 



(30) 



Substitute a\a and a!^ in eqs. ((27)) and ((29|) into ea.([23)). 
the amplitude A\ can be reduced to 



M = —^frjc / du(t){u,^) dx dymi (31) 

167r^ Jo Ja Jo 

'2 I + {2u - l)x - y 1 + (1 - 2u)x - y' 

ai 



02 



03 



• '"(fci)75w(fc2) 

, l^x-y-z 

dz 

04 



2 u{ki)a'^''v{k2)ep.i,pak^k2 



The contribution of Fig[2I^b) can be calculated in the 
same way. It can be finally shown that the contribu- 
tion of Fig[2Ub) A2 is the same as that of FiglSJ^a) if the 



condition 0(u, fi) = 4>{u, fi) maintained. This condition 
is well satisfied by the wave function of r]^. 

Then one can obtain the total contribution of the loop- 
diagrams 



with 



Ai=Ai+ A2 (32) 
= Ciu{ki)j5v{k2) + C2u{ki)a^''v{k2)e^^p„k1k2, 



Cl 



l-x 



-^^/t;,™; j du(j){u,fi) J dx J dy (33) 

2 l + (2u-l)a;- y I + {1 - 2u)x - y 
oi 



02 



03 



Q2g4 rl 

^2 = -i^^c"^i / du4){u,^i) 

dx I dy I dz- ^ 



(34) 















02 



where ai, ^ are given in eqs. ((28)) and pO|) . 

The coefficients ci and C2 are calculated analytically. 
Before performing the integrations over Feynman param- 
eters, let us make the phase convention in the integrals. 
The imaginary angles 0's for any imaginary quantities 
are restricted in the range — tt < 9 < tt. Then the loga- 
rithms appeared during the integration have a cut along 
the negative real axis. For any two imaginary quantities 
a and &. the following relations about the logarithm are 
held M 



ln(a6) = ln(a) + ln(6) + T]{a, b), 
r]{a,b) 2TTi{e{-Ima)e{-Imb)e{lmab) 
- e{Ima)e{lmb)e{-Imab), 



(35) 



where 6{x) is the unitstep function. 

The analytical results for the Feynman parameter in- 
tegration are obtained. To express the results conve- 
niently, some quantities are defined as Vc = m^/m^^^, 
n = mf/ml^, a = , and /3 = \/l -4r(. The 

result for ci is 



Cl 



-^^[h{u) + f2{u) + Uu)], (36) 



where the functions fi{u), /2(u) and fsiu) are 



/i (u) = , , \ ,^ [2re ln(r,) + u(l - 2r, - 3m + 2^^) ln(r, - {I - uf - ie) + [I - u) (37) 

(1 — u)u{rc + (1 — u)u) 

X (wln(-rc + (1 - uf + ie) - (2rc + u~ 2u^)\n{rc - u^ ~ ie) + uln(-rc + u^ + ie))], 
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/^H - {2(-l + 4n)[ln(2)ln(i^) -ln(2)ln(l±|) + 2Li.(-JL_) (38) 



- 2/31n(i^) - ln(2) ln([:j:|) + fi' ln(2) In(i^) + ln(2) h\^) " ln(2) In(i^) 

+ ^'Li2(Y^)]} + ii^(-^i) + i'l H-t'i) h H-t,) + 4 ln(-4) 

+ (-1 + ii - /3) ln(l - ti + /3) + (1 - + /?) ln(l - t[ + /3) + (-1 + t2 - /3) ln(l - + /3) 



+ (l-t'2+/3)ln(l-t^ +/?)] + 



(l-4r;)2^ Q;(ii+Q;) Q;(i'i + a) Q;(i2 + a) ^(ij, + a) 
ln(a) ln(a) ln(a) ln(a) (1 - h + /?) ln(l - h + (3) (1 - + /3) ln(l - i'^ + /3) 
ii+a i'l+a ^2 + a t'^ + a {ti + a){l + a + P) {t[ + a){l + a + (3 

(l-t2 + /?) ln(l-f2+/3) (1 - t'2 + /3) lii(l -t'2+f3) ln(l + a + /3) _ ln(l + a + /3) 
(i2 + a)(l + a + /3) (t^ + a)(l + a + /3) <i + a i'l + Q' 



ln(l + a + /3) ln(l + a + /3) ^ , 1 , 2r; + (-1 + VT^4f[)u ^^^,^ , 

+ ^^T^ ^T^l + /33(2.,+/3 + /32) ^"( ~2n + (l + Vl-4n)u ^[^(^ + 

+ /? + + 4n(2n + /? + ln(2n) - 4n(2n + /3 + /S^) ln(2n + /? + /J^)] 

- 7r4^H-ii)M-^) -ln(-t;)ln(-^) +lnH2)ln(-^) 
H-t'2) H-r^) - ln(l -h+l3) ln(i±^) + ln(l - t\ + (3) ln( ^ + " + ^ ) 

- ln(l - t2 + /-i) ln(i±^) + ln(l - t^p + /?) ln(i±^) + U^i-^) U^i-^) 

12 + a t2 + a ti + a l\ + a 

+ Li2(-— — - L12 -7-— - L12 — — + L12 — -p— ) - Ll2( — — 

12 + a 12 + a ti + a ti + a 12 + a 

, + ti, — *7r 1 

+ ^^^^ + . + ^ + 2n(n + (-i + .H [^' ^'^(-^'^ 

+ (r, + (-1 + u)u) ln((l - it)u) - ri ln((-l + - (r, - u + u^) In(w^)] | , 



r 



where 



ti = 



2n(Vi 


-4ri 4 


-u) 4- 


2n(l- 


^) 


2ri 


-(1 + 


VI- 


- 4r/)u 




2n(Vi 


- 4ri 4 


-u)- 


2n(l- 


^) 


2r, 


-(1 + 


VI- 


^ 4r/)u 




2n(Vi 


-iri - 




2n(l - 




2r,- 


h(-l4-^/^ 


- 4n)ii 




2n(Vi - 


- 4n - 


- u) + 


2n(l- 




2n - 


K-i- 




- 4n)ii 





+ ie, (39) 
(40) 
+ ie, (41) 
- ie. (42) 



The function Li2(z) is the polylogarithni fmiction, which 
is defined as Li2(z) = (jt '"*-"^^"*-* . Finally, the third 
function is 

h{u)^f2{l-u). (43) 
The analytical result of the 3-fold Feynman parameter 



integration in C2 (eq.([34|)) is tedious, it is not presented 
here. Some formulas used in the procedure of integration 
are given in Appendix B. 

It is not difficult to calculate the contribution of weak 
interaction, where the virtual Z'^ acts as the intermedi- 
ate state (Figl^fc)). The contribution of FigHJc) to the 
amplitude of rjc —i- is 



-cVcczclvc), (44) 



— m 

where p = fci + ^2, and Vji^ and VccZ are the coupling 
vertices of IIZ and ccZ, respectively. 



Vnz = 



-4 sin^ 



w 



1 



4 sin 9w cos 6w 4 sin 9w cos 9w 



75 



V-c 



ccZ 



8 sin^ 



1 



12sin6'vt'Cos6'n/ 4sin0vFCOs6'vi/ 



75 
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with 9w being the Weinberg angle. 

Contracting the leptonic final state with the leptonic 
fields operator and using the definition of the decay con- 
stant of rjc, one can finally obtain the amplitude 



Next I will go on to the numerical discussion on rjc 



(4 sin 9w cos Ow)^ — 



•^"^ u(fci)75i^(fc2). (45) 



Then the amplitude of the leptonic decay r]c ^ l^l is 
the sum of the loop and tree diagrams 

Ai7],^l+l-)=Ai+Azo. (46) 

The decay width of this decay is 

r(r?e ^ l+n = ^\MVc - I'-nf^-, (47) 

'/c 

where k is the three-momentum of one of the leptons in 
the rest frame of rjc ■ The branching ratio of the decay is 
defined by 



Br{r], ^ l+r) 



Ttot 



(48) 



where F^q^ are the total decay width of the rjc meson. 

In the numerical calculation, the wave function of ric 
is taken to be [131 



(u, fj, ~ TOc) = N 4u(l — u)exp 



f3 



4u(l - u) 



(49) 



where N is the normalization factor, the parameter /3 = 
3.8 ± 0.7. The value of the decay constant of rjc is f,j^ = 
0.346±0.017GeV[i3,the mass of c quark = 1.4 GeV, 
the total decay width of 77c is r^t = 26.7 ± 3.0 MeV [l|. 

With the parameter inputs and the wave function of 
r]c given above, the prediction to the branching ratio of 
Vc 77 is 



• 77) = 2.43_o 34 



X 10 



-4 



(50) 



where the uncertainty comes from the uncertainty of the 
total decay width of rjc, the uncertainty of the decay con- 
stant fri^ and the uncertainty of the parameter /? in the 
wave function. The main contribution comes from the 
uncertainties of the total decay width and decay constant 
of 7]c, which are about 10% of the central value, while the 
error caused by the parameter /? is small, which is only 
about 1%. 

The prediction in ea. (j50p is in good agreement with 
the experimental data 



S°"P(r?c->77) 



2.41J] X 10"". 



l+l- 



The loop integral is infrared divergent if the lepton 
mass is zero. This can be shown by taking the asymptotic 
limit of r; ^ for the analytical results of loop-integral- 
functions fi{u), f2{u) and fsiu) in eqs. ([37]) . ((38|) and 
(j43|l . The most singular term behaves like ;i ln(r;). Con- 
sidering the decay amplitude in eq. (|3ip is proportional 
to mi ~ due to the helicity suppression, the most sin- 
gular term in the decay amplitude behaves like ln(r/). 
The curves for the coefficient ci as the lepton mass ap- 
proaching zero are shown in FiglHl Only when the lepton 
mass mi < 0.2me, here me is the mass of electron, does 
the numerical result begin to be severely affected by the 
infrared divergence. As the lepton mass not less than 
0.2me, the numerical value is reliable. The situation for 
C2 is similar to ci. Therefore, for the cases of mi = me 
and mi = (m^ the muon mass), the loop integral is 
not affected by the infrared divergence. 

The values of ci and C2 are presented in Table HI 



(xlO-^) 

1.75- 

1.5 
1.25 
1 

0.75 
0.5 
0.25 



Re ci 




(xlO-^) 

2p 

1.75- 

1.5 
1.25 

1 

0.75 
0.5 
0.25 



0.2 0.4 0.6 0.8 1 

, , mi /me 

(a) 



Im Ci 



0.2 0.4 0.6 0.8 1 
mi/me 



(b) 



FIG. 3: Behavior of the coefficient ci as mt 0. (a) the real 
part of ci; (b) the imaginary part of ci. The horizontal axis 
is for the ratio of mi/nie. 



(51) 
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TABLE L The coefficient ci and C2 in tfie decay ampiitudes of 77c e^e and 77c — > M^A* 





Cl 


C2 


77c e+e 


5.90 X 10^'' + 1.12 X lO^'i 


-2.75 X lO"'' - 1.43 X lO^^'i 


77c ^ /^"""m" 


3.20 X 10"" + 1.02 X lO^^'i 


-5.53 X 10" ' +5.36 x IQ-'^i 



The branching ratios are given in Table HIl The lep- 
tonic decays of rjc are dominated by the two-photon loop 
diagrams, the contribution of tree diagram modifies 
the decay rates by less than 1%. 



TABLE IL The branching ratios of rjc —> e^e" and 77c —>■ 
fi'^ . The column "Only loop" means the contribution of 
only loop diagram, without the contribution of tree di- 
agram, while "Loop-t-^"" denotes the contribution of both 
loop and tree diagrams. 





Only loop 


Loop-fZ" 


77c -> e+e 


4.77t;!-^ X 10-^^ 


4-74lo;e6 X 10"'^ 




6.41±i-;« X 10-« 


6-391:^;^;; X IQ-'" 




(a) (b) 

FIG. 4: Diagrams for 77c in physics beyond the stan- 

dard model, (a) leptoquark contribution, and (b) light pseu- 
doscalar Higgs contribution. 



The branching ratios of 77c — > are very tiny within 
the frame work of the standard model. Charm physics in- 
cluding leptonic, semileptonic and hadronic charm decays 
wiU be studied at BESIII About 10 x 10^ J/tA 

events will be accumulated with one year designed lumi- 
nosity. Considering the branching ratio of J/t/i ^rjc 
being (1.3 ± 0.4)% Q, about 1.3 X 10 rjc mesons can 
be produced through the radiative decays oi J/il) jrjc. 
Assuming the detection efficiency is about 20%, then the 
sensitivity of the measurement of 77c — ^ can reach 

3.8 X 10~*. The standard model prediction to the lep- 
tonic decay of 77c is just below the sensitivity of 
BESIII. However, the low decay rate within the standard 
model makes the leptonic decay of 77^ l^l^ sensitive to 
physics beyond the standard model, such as the existence 
of leptoquark bosons. Fig. Ufa), or the light pseudoscalar 
Higgs boson, Fig. HJ^b) . Especially when the mass of the 
light pseudoscalar Higgs lies near the mass of 77c meson. 
Fig. m^b) can significantly enhance the leptonic decay 
rate of 77c. Measurement of 77c decay, especially 
for 77c fJ''^ fJ-~ channel, can give information of physics 
beyond the standard model. 

In summary, the radiative and leptonic decays of 77c 
77 and 77c have been studied consistently within 

the framework of the standard model. The theoretical 
prediction of the decay rate of 77c — > 77 is in well agree- 
ment with the experimental measurement. For the lep- 
tonic decays of 77c — > ^J■^^J'^ ^^nd e^e~, within the frame- 
work of the standard model, the loop diagrams of elec- 
tromagnetic interaction dominate, the tree diagram of 
weak interaction involving z'^ propagator can only mod- 
ifies the decay rates by less than 1%. The decay rates of 



77c /^^M^ and e+e~ arc tiny, which makes them sen- 
sitive to some new physics beyond the standard model. 
Measurement of these decay rates may shed light on the 
existence of new physics. 
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Appendix A 



The wave function of a meson M composed of a quark- 
antiquark pair 51^2 can be defined through the ma- 
trix element {0\qi{x)aq2iy)f3\M) , where summation over 
color degrees of freedom and a gauge invariant factor 
Pexp{i J dz^T"- A'^{z)} between qi{x) and (72(2/) are in- 
dicated. 

Using the equality of Fictz transformation 



SikSij = -^SijSik + ^(75)^(75);/: (52) 

+ \h^^)^jh^)lk - ^(7A<75)y(7^75);fc 

+ ^(o-^..75)«j(o-^''75)ifc, 
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the matrix element can be written as 

mi{x)c.q2{v)(3\M) 

= {^\qi{x)pq2{y)a\M)5pa8i3„ 



(53) 



= {^\(ll{x)pq2{v)a\M){^5p„5pa + ^(75)p<T(75)/3a 
1 1 



contributions of the wave functions of twist-3 are sup- 
pressed due to their spin structure. Therefore they are 
safely dropped. 

Appendix B 

Some formulas used in the calculation of the 3-fold Feyn- 
man parameter integration are 



+ ^{<yii.ylb)pa{<y^''lb)pa) 



dx 



1 



1 



\qi{x)lBq2{y)\M){-fr,)[ja 



ax^ -|- 6a; -|- c 
- \J\P- ~ iac ] 



ln(2aa; + b (58) 



{0\qi{xhpq2iy)\M){r)pa 
(0|gi(a;)7M75g2(y)|M)(7^75)^„ 

(0|gi(a;)a^.7592(y)|A/)(a^''75)/3a. 



da; 



ln(2aa; + b+ \/p~4ac 
bx + 2c 



{ax^ +bx + c)^ 
b 



dx 



For a pseudoscalar meson A/, the matrix elements of the 
scalar and vector current are all zero. For the nonzero 
matrix element the wave function can be defined as po| 

{0\qi{xh^l5q2{y)\M) ^ tfMPp due-^^^^-^+^f </)(«, A^), 

Jo 

(54) 



62 - 4ac 
For the integration of the type 

niy + n2 



(6^ — iac){ax'^ 
1 

ax^ + 



bx 



(59) 



bx 



dy-, 



(0|5i(x)z75g2(2/)|M) 



M 



mi + 1712 Jo 



(55) 



(aiy2 + biy + Ci){a2y^ + 622/ + C2) '' 

the integrand can be decomposed as 

niy + 712 

(ai?/2 + biy + Ci)(a2y2 + 622/ + C2) 
hy+h , hy+h 



(60) 



(61) 



{0\qi{x)ap^-f5q2{y)\M) = i 



fAim 



M 



iPp.Zu - PuZp.) (56) 



mi + m2 
6 ' 



aij/^ + biy + ci 02?/^ + 622/ + C2 

where the coefficients ki^2 and ^1^2 can be solved within 
the above equation. Then the integration of eq. (j60p can 
be reduced into a simpler form. 
For the integral like 



where z ~ y ~ x, fi is the energy scalar where the pcr- 
turbative and non-perturbative dynamics of QCD can be 
factorized, /m the decay constant, pp the momentum of 
the meson M, and ttim, "t-i and m2 are the masses of 
the meson M, quarks qi and q2, respectively. The wave 
function 4){u^^) is of twist-2, and 4>p{u, ^) and 0cr(u,/i) 
are of twist-3. 

Then the matrix clement {0\qi{x)aq2{y)t3\M) becomes 



dx- 



kx + I 



1 



:ln(.g.a; + 6± ^62; + /), (62) 



ax^ + 6a; + c -^ea; + / 

the integration can be performed by making the variable 
transformation t = y/ex + f. While for 

kx + l 1 



dx 



aix^ + bix + ci ^a2X^ + 62a; -I- C2 



(63) 



ln(5a; -f 6 ± ^ a2x'^ + 62a; -I- C2), 



]qiix)aq2{y)t3\M) 

--Jm I due-''^^P--+''P-y\i>-f,ct>{u,ii) 



(57) 



with 02 > 0, the variable transformation shall be 



020;^ 



'A/ 



mi H- m2 



75(^p(u,^) 



1 m 



M 



2 mi -|- m2 



{Pti-Zu 



62a; -I- C2 = t — y/a^x. 



^^ln(t_a) ^ 



t - a 



-PuZp)apy 75j/3a. 



-Li2 - 



a — b 
t - a 



(64) 



a 



For the decay processes of rjc considered in this work, 
only the leading-twist wave function is considered. The 



where Li2(2:) is the polylogarithm function, which is de- 
fined as Li2(z) = 
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